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,  ABSTRACT  AND  SUMMARY 

This  report  presents  a  mechanism  for  calculation  of  the  optical  Index 
structure  constant  from  bulk  meteorological  measurements  (air  and 
water  temperatures,  wind  speed,  air  humidity)  made  over  the  ocean.  From  this 
and  a  formulation  of  the  near  field  optical  propagation  from  an  Incoherent 
point  source,  the  expected  line  spread  function  for  near  the  sea  surface  lines 
of  sight  Is  calculated.  Finally,  three  series  of  photographic  line  width 
observations,  made  during  unstable,  neutral,  and  stable  atmospheric  buoyancy 
conditions,  are  compared  with  line  widths  expected  from  the  calculations. 

Mean  values  of  observed  results  agree  with  the  calculations  within  about  25 
>ercwt"for  conditions  of  strong  atmospheric  instability.  Since  comparisons 
were  made  over  a  variety  of  target  ranges,  wind  speeds,  and  camera  elevations, 
it  is  concluded  that  the  formulations  derived  are  accurate  to  at  least  £  25  ^ 
percent  for  conditions  of  strong  atmospheric  Instability. 

The  observations  made  under  conditions  of  neutral  and  stable  atmospheric 
buoyancy  were  Instrumental 1y  limited  to  an  angular  subtense  of  2  to  4  arc 
seconds.  These  measured  line  width  values  fall  between  1.9  and  0.6  times  the 
values  computed  from  meteorological  measurements. 

The  variability  of  image  spreading  observed  across  the  target  fields 
points  out  the  difficulty  of  employing  a  single  mechanism  for  correction  of 
atmospherically  induced  image  smear. 
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I.  INTRODUCTION  AND  BACKGROUND 

The  information  collection  capability  of  photography  with  a  camera  located 
near  the  sea  surface  is  limited  by  the  luminance  characteristics  of  the 
target,  by  the  characteristics  of  the  atmospheric  path,  and  by  the 
characteristics  of  the  camera  in  its  local  environment.  Significant  losses  of 
information  arise  from  the  loss  of  contrast  and  picture  smearing  resulting 
from  the  characteristics  of  the  atmosphere  along  the  line-of-sight. 
Relationships  predicting  the  loss  of  contrast  are  relatively  well  known;! 
however,  the  relationships  between  picture  smearing  (attenuation  of  optical 
modulation  transfer),  sea  state,  height  of  the  line-of-sight  above  the 
surface,  and  meteorological  parameters  have  only  recently  been  developed.  In 
particular,  the  formulation  of  the  relationship  of  the  strength  of  the 
turbulence  Cfl  to  gross  meteorological  parameters  in  the  near-the- 
ocean-surface  environment  has  been  under  development  only  since  1974,2  an(j 
no  measurements  of  Cjj  have  been  made  closer  than  one  meter  from  the  sea 
surface.  This  independent  research  effort  is  one  of  the  first  which  relates 
meteorological  variables,  turbulence,  and  picture  smearing  for  a 
near-the-ocean  line  of  sight.  The  ultimate  goal  is  to  be  able  to  predict  the 
smearing  due  to  turbulence,  given  details  on  the  camera  and  target  geometries 
and  the  gross  meteorology  at  a  particular  site. 

In  a  previous  NUSC  effort-*  the  recording  camera  was  located  in  a  tent  on 
a  beach.  The  target  was  mounted  on  a  ship  stationed  from  2  to  3.7  km.  away. 
The  effective  turbulence  level  for  the  integrated  line  of  sight  was  much 
larger  than  the  turbulence  level  measured  strictly  over  the  water.  The  result 
was  attributed  to  turbulence  generated  by  breakers  and  by  the  shore. 

In  the  current  effort,  the  camera  equipment  was  located  on  an  offshore 
tower  to  provide  a  full  over-water  line  of  sight,  and  meteorological  measuring 
equipment  was  located  nearby  (c.f.  Fig.  1).  The  meteorological  measurements 
performed  for  this  study  include  primarily  macro  (bulk)  meteorological 
measurements.  A  few  micro  meteorological  measurements  were  obtained  for  use 
in  comparison  calculations  of  c£. 

In  this  memorandum  a  mechanism  for  deriving  values  of  the  index  structure 
constant,  Cfj,  from  values  of  bulk  meteorological  parameters  is 
presented.*  From  this,  using  near  field  optical  propagation  formulations, 
expected  values  of  the  width  of  the  optical  line  spread  function  for 
Incoherent  imaging  are  obtained.  These  are  compared  with  averages  of  measured 
line  spread  function  widths.  Criticisms  of  experimental  results  and 
conclusions  follow. 

The  meteorological  portion  of  the  study  was  guided  by  personnel  of  the 
Meteorology  and  Physics  Department  of  the  Naval  Postgraduate  School  of 
Monterey,  California.  Measurements  of  variables,  data  reduction,  and  optical 
propagation  portions  of  the  study  were  carried  out  by  NUSC  personnel. 

II.  TURBULENCE  THEORY 
a.  Structure  Functions 

In  the  nature  of  turbulence  produced  by  shear  flow,  the  initial  power 
input  is  into  large  eddies.  This  power  is  then  transferred  from  large  to 


*  Details  of  the  approach  have  been  given  by  Davidson,  et  alia,  in  Reference  4. 
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smaller  eddies  with  little  energy  loss  until  It  Is  finally  converted  Into  heat 
by  viscous  forces  when  the  eddy  sizes  approach  the  Inner  scale  of  turbulence 
lo  (with  dimensions  of  a  few  mm.).  When  the  shear  Is  produced  by  friction 
In  flow  parallel  to  a  boundary,  the  scale  size  L0  of  the  input  eddies  at  a 
distance  z  from  the  boundary  Is  approximately  equal  to  z.  For  distances  r 
lying  in  this  Inertial  range,  i.e.,  L0  >  r  >  10»  we  may  define  a  velocity 
structure  function  Dv  giving  the  mean  square  velocity  Y  difference  for  two 
points  separated  by  a  distance  r, 

DV(M>  *  <  CV(0)  -  V(r)]2  >,  (1) 


where  for  simplicity  the  definition  is  expressed  in  scalar  fashion,*  and  where 
over  distances  of  the  order  of  r,  the  turbulence  is  Isotropic.  For  distances 
r  within  the  inertial  range,  the  Kolmogoroff  or  Von  Karman  spectrum  of 
turbulence  element  sizes  leads  to  the  relationship: 


Vr) 


(2a) 


It  has  been  shown  that  when  gradients  of  passive  additive  elements,  e.g., 
temperature  e  or  humidity  Q,  exist  in  the  shear  flow,  the  spectra  of  element 
sizes  lead  to  similar  relationships,  i.e., 

D0(r)  =  .  r2/3  (2b) 

and 

dq( r)  =  Cq  •  r2/3  (2c) 


respectively,  where  the  proportionality  constants  C  are  "constant"  over 
regions  having  dimensions  of  the  order  of  fractions  of  meters,  but  may  vary 
slowly  over  larger  distances.  Since  the  changes  In  optical  Index  n  are 
related  to  air  density  and  humidity  changes,  we  have  a  similar  relationship 
for  this  variable 


Dn(r)  *  Cn 
n  n 


(3) 


Cfl  has  been  called  the  optical  Index  "structure  parameter." 


*  Properly  used,  Dv  is  a  second  rank  tensor. 
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b.  General  Macrometeorologlcal  Relationships 

For  calculations  of  the  strength  of  the  optical  structure  parameter  C2 
from  bulk  meteorological  parameter  measurements,  the  following  input  variables 
are  used: 


T(z) 

St,i 


height  above  ocean  surface  (m. ) 

air  temperature  (OK.) 

water  surface  temperature  (°K.) 

humidity  mixing  ratio  (gms.  water/kg.  air) 

humidity  mixing  ratio  at  sea  surface  (100%  r. 

wind  speed  (m./sec.) 

atmospheric  pressure  (millibars) 


h.  at  T  =  Tn) 


From  these  one  uses  the  functions: 

9  (z)  *  potential  temperature  =  T(z)  +  .Olz 
A  (2)  *  (Z)  -  Tq 

Aq( z)  s  q(z)  -  q0 
£  »  z/L 


where 

L  =  Monin-Obukhov  scaling  length  (m. ) 
is  a  constant  obtained  from  formulations  of  a6,  aQ,  and  U. 

Since  the  formulations  to  be  used  differ  markedly  with  the  convective 
stability  of  the  atmosphere,  we  proceed  to  define  stability  by 

A0(z)  +  6.1  x  10-4T(z)  Aq(z)  >  0  :  stable  atmosphere 

a 0( z)  +  6.1  x  10_4T(z)  Aq(z)  <  0  :  unstable  astmosphere 

A0(z)  +  6.1  x  10-4T(z)  Aq(z)  =  0  :  neutral  atmosphere 

Following  the  approach  of  Ref.  4,  It  is  then  possible  to  utilize  values  of 
T,  T0,  q,  U  and  P  measured  at  a  height  of  convenience  (e.g.  5  to  20  m. )  to 
obtain  a  continous  formulation  giving  C fj(z).  The  necessary  steps  are 
outlined  In  Appendix  I.  The  resulting  expression  is 

C2  =  6.2x10'9.P2.T"4-[t|+0.070T*Q*  +  0.0019  qJJ  z"2/3  f(z/L),  (4) 


where  P  is  pressure  In  mb 

T  Is  absolute  temperature 

T*,  Q*  and  L  are  computed  in  accordance  with  Appendix  I 
f ( z/L >  «  4.9  [1+2.4  (z/L)2/3]  for  a  stable  atmosphere 
and 

f ( z/L )  •  4.9  [1-7  z/L]-2/3  for  an  unstable  atmosphere 


(5. a) 


( 5.  b) 
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With  this  formulation,  we  proceed  to  determining  the  effect  of  Cft  upon 
imaging  over  a  slant  path  above  the  sea  surface. 

III.  OPTICAL  PROPAGATION  AND  IMAGING  THEORY 


a.  Effective  Refractive  Irregularity  Strength  Over  the  Line  of  Sight 

To  treat  the  propagation  of  imaging  light  from  the  target  to  the  camera  in 
the  optical  "near  field,"*  we  assume  that  the  target  consists  of  a  series  of 
incoherent  point  sources  of  spatially  varying  Intensity.  Imaging  light  from 
each  source  then  acts  independently  and  travels  through  the  turbulent  medium 
to  a  camera  at  a  range  R.  With  a  point  source,  turbulence  elements  affect  the 
image  in  accordance  with  the  5/3  power  of  the  distance  from  the  source. 4, 5 
Thus,  the  effective  product  of  the  refractive  irregularity  strength  and  the 
path  length  becomes: 


C^(eff)-R 


(x) - (x/R) 5/3*dx , 


(6) 


where  x  is  the  distance  from  the  source.  To  simplify  our  analysis,  we  assume 
that  the  turbulence  is  uniform  over  the  range,  i.e.,  the  turbulence  is 
statistically  uniform  at  ar\y  particular  height  h  above  the  water,  but  the 
average  value  varies  with  h,  or 


C^(eff ) *R  =  f  [h(x)]*(x/R)5/3*dx.  (7) 

0J 


In  terms  of  the  geometry  of  Figure  1, 


h  «  ht  -  x(ht  -  hc )/R  -  ht(l  -  ax/R), 


(8) 


where 


ht  *  height  of  target 
he  »  height  of  camera 
a  *  1  -  hc/ht* 


For  the  case  of  a  stable  atmosphere,  if  we  use  Equations  (4)  and  (5. a)  the 
height  dependence  of  Cf|  may  be  expressed  by 


elm 


C2(h  ) 
n  m 


(9) 


*  In  the  "near  field"  of  optical  propagation  wavefront  distortions  of  the 
receiver  are  primarily  optical  phase  distortions;  whereas  In  the  "far  field" 
optical  amplitude  distortions,  leading  to  image  scintillation,  also  occur. 
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where  hm  Is  the  height  at  which  Cfj  is  measured  or  derived  from  bulk 
meteorological  measurements.  Substituting  Cjj  from  Equation  (9)  in 
Equation  (7)  and  letting  x/R  =  u  gives: 


2  R  ‘  Cn(h J  ‘  h mZ/3  TM 

C2(eff)  •  R  3  n  m  m 
n 


[1  +  2.40L/L)*7*]  L 


m 


7 u5/3" +  U* 7  TT^] 


(10) 


=  R ;  ^2/3cX) 
[ 


- oi-  .  1"^-  +  -i-  •  I  ( a,  oil 

1  ♦  2.4(hm/L)2/3]  L2/3  ht2/3  j 

m  z 


where  the  function: 


I(a,B)  =  (1  +  7b) 


2/3 


/ 


u5/3du 


{(1  -  au)  [1  +  7e(l  -  au)]}' 


m 


(ID 


wi  th 


a  =  1  -  hc/ht 

B  *  -ht/L, 

is  given  by  the  graph  of  Figure  2. 

For  the  case  of  an  unstable  atmosphere.  Equations  (4)  and  (5.b)  give  the 
height  dependence  of  Cfl  as: 


C2(h)  -  C2(hm 
n  n  ro 


h id  -  7h m/U"l2/3 
m _ m 

h(l  -  7h/L)  J 


(12) 


Substituting  for  h  from  Equation  (8)  in  Equation  (12)  and  also  Equation  (12) 
in  Equation  (7)  gives  for  Cjj(eff)*R 


C2(eff  )*R 


R*W' 


2/3 


I(o,b)  . 


(13) 
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Values  of  I (a,  e)  may  be  obtained  similarly  from  the  curves  of  Figure  2. 

b.  Mutual  Coherence  Function 

Yura7  gives  an  expression  for  the  averaged  mutual  coherence  function 
(MCF)  of  light  from  an  incoherent  point  source  propagated  through  a  turbulent 
atmosphere  to  an  aperture  of  diameter  D  as: 

Mst(p,R)  =  exp^  (p/p0)5/3  -  $  •  ( p/D)1/3j  j>-  (14. a) 

where 

Mst  =  the  averaged  short  term  (short  exposure)  MCF 
p  =  xfF  =  wavelength-camera  focal  length-spatial  frequency  in 
image  plane 

Po  =  long  term  (long  exposure)  coherence  distance 

_  2.91  r2,  wV  _  -3/5 
-  ~2  Cn(eff )‘^— y  -R 

(14. b) 

=  1.60  x  10"14C2(eff)-R  '3/5 

n 

x  =  6  x  10“7  m.  (approx. ) 

«  =  multiplier  in  the  range  from  0.5  to  1.0  depending  upon  the  target 

range.  0.5  corresponds  to  far  field  imaging;  1.0  corresponds  to 
near  field  imaging.*  In  our  case,  we  used  the  multiplier  0.75  in 
all  the  calculations,  except  in  the  1980  series-  Since  the 
atmosphere  was  neutral  then,  imaging  was  near  field,  and  6  =  1.0. 

c.  Camera  Line  Spread  Function 

Equation  (14. a)  is  also  in  the  form  of  a  statistically-averaged 
modulation  transfer  function  MTF  for  the  cone  of  atmosphere  from  the  receiving 
aperture  to  a  target  point.  To  obtain  the  MTF  for  the  full  recording  system 
MTFSyS,  we  multiplied  the  atmospheric  response  by  the  (approximated) 
frequency  response  H  of  an  ideal  telescope  having  an  86  mm  window  and  a  25  mm 
obscuration  and  by  an  approximation  to  the  frequency  response  J  of  Pan-X  film 
obtained  from  a  handbook  graph.  Then  substituting  the  value  6  x  10-7m  for  x 
and  1.47  m  for  f  gives  for  the  overall  response 

MTF$ys(F,po)  =  exp[-8.8  x  10 "7 - ( F/pQ ) 5/3 • ( 1  -  0.0219  6.F1/3)]-H(F)-J(F),  (15) 

2 

where  pQ(Cn,R)  was  defined  in  Equation  (14. b). 


pq  <<  xR.  imaging  is  far  field.  If  p0  >>  xR,  imaging  is  near  field. 
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H  =  F/ 60, 000 
H  =  1.0/1. 4  -  F/ 140,000 
H  5  0 
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If  0  <  F  «  30,000 
if  30,000  <  F  «  100,000 
if  F  >  100,000 


and 

J  =  [l  +  3(F/53,000)2]‘1/2, 
where  the  units  of  F  are  cycles/m. 

To  compare  empirical  with  theoretical  results,  one  might  Fourier  transform 
the  observed  image  spread  functions  and  compare  the  averaged  transforms  with 
Equation  (15).  Rather  than  perform  the  multiplied  measurements  and  transforms 
which  this  would  require,  we  Fourier  transformed  Equation  (15)  to  obtain 
expected  average  line  spread  functions  for  a  variety  of  p0  values  and  for  6 
=  0.75  and  1.0.  Figure  3  gives  smoothed  curves  showing  the  values  of  line 
widths  (1/e  times  peak  amplitude)  obtained  from  these  transforms.  These  were 
used  for  determining  expected  line  widths. 

IV.  EXPERIMENTAL  ARRANGEMENT 

The  line  of  sight  from  the  target  to  the  camera  was  as  shown  in  Figure 
1.  The  camera,  a  “Questar"  telescope  coupled  to  a  Nikon  35  mm.  camera  back, 
was  attached  to  a  tower  corner  post  at  a  variety  of  heights  1^-  above  mean 
water  level.  The  target  consisted  of  a  series  of  five  1000  watt  quartz-iodide 
bulbs,  each  having  a  luminous  area  of  160  by  1.0  mm.  The  target  was  placed  at 
a  variety  of  heights  h^  at  ranges  R  from  460  to  5000  m. 

The  meteorological  measuring  equipment  was  located  on  an  adjacent  corner 
of  the  tower  in  the  1978  observation  series  and  on  a  mast  rigged  on  the  bow  of 
the  freighter  YFRT-287  in  the  1980  and  1981  observation  series.  The  freighter 
was  anchored  460  m  from  the  tower,  just  downwind  from  the  camera's  line  of 
sight.  The  equipment  consisted  of  an  anemometer,  wind  vane,  dew  point 
apparatus,  air-water  temperature  difference  measuring  apparatus,  and  water 
temperature  measuring  apparatus.  The  in-air  sensors  were  located  5.7  m.  above 
the  water  level  in  the  1978  series  and  7.0  and  8.1  m.  above  the  water  level  in 
1980  and  1981.  Water  temperature  measurements  were  of  the  "bucket"  variety. 

In  the  1980  observation  series,  a  micrometeorological  (differential  air 
temperature  measuring)  apparatus*  and  recorder  were  mounted  on  the  Brenton 
Reef  Tower.  By  means  of  differential  temperature  measurements  made  by  this 
instrument  in  the  air  stream,  values  for  the  temperature  structure  function  of 
the  air  could  be  obtained  directly.  And  these  were  simply  convertible  to  the 
optical  index  structure  function  which  was  needed  for  describing  and 
calculating  optical  image  propagation.  This  apparatus  was  used  on  the  data 
taking  effort  of  16  May,  1980.  The  results  obtained  from  this  apparatus  were 
compared  with  values  of  Cfj  deduced  from  simultaneously-made  bulk 
meteorological  measurements.  When  we  attempted  to  use  the  micrometeorological 
measuring  apparatus  on  the  primary  data  taking  exercise  of  6  June,  1980,  the 
fragility  of  the  probes  and  the  lack  of  a  sufficient  number  of  spares  limited 
our  data  taking  to  four  single  measurements. 


ontel ,  Inc.,  Model  MT-2  microthermal  probe 


LONG  TERM  COHERENCE  DISTANCE 
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Unfortunately,  the  air-water  differential  temperature  apparatus  and  the 
mlcrometeorologlcal  apparatus  failed  to  operate  In  the  1981  data  taking 
series.  The  air-water  differential  temperature  apparatus  was  replaced  with  a 
thermometer  located  on  the  shady  side  of  the  ship's  mast  and  a  thermometer 
kept  in  a  bucket  in  the  water. 

Accuracy  limitations  of  the  sensors  were  approximately  as  follows: 

Static  camera  resolution  (with  Pan-X  film)  .  2.5  arc  sec.  (equlv.  to 

18  microns  on  the  film) 


Image  smearing  due  to  target  motion . 0.04  arc  sec.,  max.,  in 

1/250  sec.  exposure 


Luminous  target  geometrical  line  width  .  0.45  arc  sec.,  max. 

Wind  speed  .  0.5  m/sec.  starting 

speed;  running  accuracy 
+0.5  m/sec. 


Dew  point  temperature  .  +  1°C 

Air-water  differential  temperature  .  +  0.1“C  (est)  for  1978 

(air  sensor  is  aspirated  and  and  1980  series 

double-shielded  in  1978  and  1980  series)  +  0.5°C  (est)  for  1981 

series 


Water  temperature  measuring  apparatus . +.  1°C 

Differential  air  temperature . +  10-17  Cf}  units 

measuring  apparatus  “ 
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V.  DATA  TAKING  AND  PROCESSING 
a.  General 

The  three  series  of  observations  were  made  at  Brenton  Reef  Tower,  located 
just  off  the  entrance  to  Narragansett  Bay,  R.I.  The  tower's  location  vs.  the 
nearest  land  Is  shown  on  the  map  of  Figure  4.  We  attempted  to  position  the 
target  carrying  ship  In  a  direction  such  that  the  camera  line  of  sight  was 
normal  to  the  wind  direction.  Camera-target  heights  and  ranges  are  given  in 
Table  I.  In  nearly  all  cases  the  camera  exposure  time  was  1/1000  sec.  To 
ascertain  that  critically-focussed  images  were  obtained,  we  first  determined 
the  relationship  between  the  camera  viewer  best  focus  and  the  film  best  focus 
positions  in  the  laboratory.  Then  in  the  field,  at  every  target  range 
recordings  were  made  for  a  series  of  focus  positions  bracketing  the 
laboratory-determined  best  focus  position. 

TABLE  I 

GEOMETRY  FOR  FILM  RECORDING 


Observation 

Series 

Day 

Time 

Target 
Ranges  (m) 

Target 
Height  (m) 

Mean  Camera 
Height  (m) 

A 

9/26/78 

10:00-11:30 

460-2740 

6.2 

2.4 

B 

9/26/78 

11:45-12:50 

460-2740-460 

6.2 

1.4 

C 

9/26/78 

14:00-14:30 

460-1830 

6.2 

1.0 

F 

9/27/78 

12:30-12:50 

2740-460-1830 

6.2 

2.3 

G 

9/27/78 

13:25-14:00 

460-2740-460 

6.2 

0.8 

H 

6/6/80 

10:55-13:05 

990-4970-520 

9.85 

0.9 

I 

6/6/80 

14:05-15:40 

950-5145-500 

9.85 

1.7 

J 

6/6/80 

16:00-17:15 

425-5030-1000 

9.85 

3.1 

K 

5/8/81 

10:05-11:18 

500-4000-3000 

7.9 

1.1 

L 

5/8/81 

11:30-12:15 

5000-500 

7.9 

0.95 

M 

5/8/81 

12:35-14:02 

500-5000-500 

7.9 

2.0 

N 

5/8/81 

14:30-15:00 

500-5000 

7.9 

3.3 

Bulk  meteorological  data  were  taken  and  recorded  continuously  during  the 
three  expeditions.  Recorded  results  were  first  averaged  for  fifteen  or  thirty 
minute  periods;  then  the  averages  were  interpolated  to  give  data  for  the 
camera  image  recording  times. 

In  the  laboratory  the  ends  of  each  exposed  film  strip  were  given  a  step 
wedge  exposure,  and  the  film  was  processed  continuously  in  a  Versamat 
processor.  Then  by  use  of  a  microscope  and  scale,  the  image  negatives  were 
scanned  to  determine  approximate  image  widths  and  thus  to  determine  the 
best-focus  frames.  These  were  then  scanned  with  a  Joyce-Loebl  Model  3C  micro¬ 
densitometer.*  On  the  selected  frames  we  attempted  to  choose  nominally 
average  portions  of  the  five  source  images  for  scanning.  This  same  principle 


T 
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the  formulation  [(measured  width)*  -  (scanning  slit  width)*]!/*  was 

used  as  a  correction  factor  for  the  image  broadening  caused  by  finite  slit 

width. 
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was  applied  when  the  Images  were  smeared  together  because  of  atmospheric 
Irregularities.  On  any  scan  the  peak  density  value  was  obtained  from  a 
(visually)  smoothed  trace,  and  an  averaged  background  density  value  was 
obtained.  From  the  step  wedge  exposures,  the  density  vs  exposure  transfer 
(H  and  D)  curve  for  Pan-X  film  was  obtained,  and  the  equivalent  peak  and 
background  exposures  Ep  and  E5  were  found.  Then  the  l/e'th  signal 
exposure  level  Ees  was  calculated  from 

Ees  =  (Ep  -  Eb)/e,  (16. a) 

and  the  corresponding  exposure  with  background  added  was 

Ee  =  Ep/e  +  Eb(l  -  1/e).  (16. b) 

By  the  use  of  the  transfer  curve  again,  the  image  density  corresponding  to 
Ee  was  determined.  Then  the  image  width  between  l/e'th  points  was 
measured.  Finally,  for  a  particular  exposure  series  and  target  range,  the  1/e 
image  widths  were  averaged  to  obtain  mean  values,  and  expected  standard 
deviations  of  the  mean  and  of  individual  images  were  computed. 

b.  1978  Observation  Series 

The  first  series  of  observations  was  made  with  a  strongly  unstable  marine 
atmosphere  on  September  26  and  27,  1978.  Thirty  minute  averages  of  the 
recorded  meteorological  data  are  shown  on  the  eight  lower  curves  of  Figures  5a 
and  5b.  On  September  26  the  sky  was  fully  clear.  The  uppermost  curves  give 
values  of  Cfj  computed  by  Naval  Postgraduate  School  personnel  for  the 
meteorological  apparatus  height,  5.7  m. 

On  this  expedition,  a  total  of  411  exposures  were  made,  and  line  widths 
were  measured  on  73  of  these.  Figures  6  and  7  are  enlarged  prints  of 
"typical"  exposures  from  Series  A.  By  use  of  a  mask  the  light  source  portions 
of  the  negatives  were  overprinted  to  show  the  structure  in  the  spread 
functions.  Since  the  exposure  times  for  the  negatives  were  chosen  to  avoid 
saturation  of  the  target  source  portion  of  the  frame,  the  ship  and  background 
are  underexposed.  The  length  of  the  target  array  is  0.8  m. 

c.  1980  Observation  Series 

On  the  expedition  of  16  May,  1980,  the  shutter  of  the  Nikon  data  camera 
failed  to  operate,  resulting  in  no  film  exposures.  Thus,  only  micro- 
meteorological  and  bulk  meteorological  measurements  obtained  on  this  date  were 
useful.  Figure  8  compares  values  of  Cfl  computed  from  bulk  meteorological 
measurements  with  values  obtained  by  means  of  the  microthermal  bridge 
apparatus. 

The  major  results  of  the  series  were  obtained  on  the  expedition  of  6  June, 
1980.  At  this  time,  about  460  exposure  (images  of  the  target)  were  made,  and 
98  of  these  were  scanned  subsequently  with  the  laboratory's  micro¬ 
densitometer.  Each  exposure  was  scanned  about  three  times  to  obtain  a  mean 
value  for  the  spread  of  the  line  image  on  the  frame.  37  of  these  frames  were 
then  selected  as  representative  of  best  focus  conditions,  and  from  these  the 
line  widths  representative  of  the  geometrical  conditions  of  the  data  taking 
were  obtained.  Figures  9  and  10  are  representative  photographs  of  targets  at 
1000  and  4970  m.  The  quality  of  the  target  image  in  Figure  10  should  be 
compared  with  the  image  of  Figure  7  obtained  in  the  presence  of  strong 
turbulence  at  only  1830  m.  range. 
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Pigure  5.b  METEOROLOGICAL  DATA  TAKEN  ON  SEPTEMBER  27,  1978 
AND  RESULTING  COMPUTED  VALUES  OF  C2 


FIGURE  10  ENLARGED  PRINT  OF  IMAGE  OBTAINED  WITH  A  WEAKLYcTURBULENT  ATMOSPHERE 

Range:  4970  m. ;  C£ :  0.72  x  lO"15 
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Values  of  the  measured  meteorological  variables  and  computed  values  of 
for  the  data  taking  period  of  June  6  are  shown  in  Figure  11.  The  gaps 
in  the  curves  represent  the  time  interval  when  the  YFRT-287  was  transiting  to 
a  new  anchorage  location.  The  move  was  necessitated  by  the  earlier  shift  in 
the  wind  direction. 

d.  1981  Observation  Series 

On  the  expedition  of  May  8,  1981  a  series  of  about  400  exposures  of  the 
target  were  made,  and  77  of  these  were  subsequently  scanned  by  the 
laboratory's  microdensitometer.  Each  exposure  was  scanned  about  three  times, 
and  32  of  the  exposures  were  selected  as  best  focussed.  From  these  the  line 
widths  representative  of  the  geometrical  conditions  of  the  data  taking  were 
obtai ned. 

Values  of  the  measured  meteorological  variables  and  computed  values  of 
are  shown  in  Figure  12.  During  the  time  interval  from  09:30  to  11:15, 
values  of  Cfj  could  not  be  computed  from  Equations  (5. a  and  5.b),  since 
£o  >  1.25,  and  no  formulation  for  the  computation  of  C^  from  bulk 
measurements  exists  for  the  highly  (buoyantly)  stable  atmosphere  indicated  by 
the  measurements. 

VI .  RESULTS  OF  OBSERVATIONS 

a.  General 

The  three  data  taking  series  fall  into  three  categories  of  atmospheric 
stability:  unstable,  neutral,  and  stable  for  the  1978,  1980,  and  1981  series, 
respectively.  Thus  the  results  are  presented  separately  by  years  on  similar 
graphs  (Figures  13  to  24).  Details  common  to  all  the  graphs  follow. 

Arrows  on  the  curves  indicate  the  time  sequence  in  which  data  were  taken. 
Line  widths  shown  represent  the  distances  between  points  on  the  film  having 
exposure  values  with  amplitudes  1/e  times  the  peak  amplitude.  Ordinate  scales 
are  linear  throughout,  except  for  the  "b"  series  (1978  only)  which  uses  a 
logarithmic  ordinate  to  illustrate  the  percentage  inaccuracy  between  computed 
and  measured  values  of  line  width.  The  abscissa  of  each  graph  is  plotted  on  a 
(range)0*6  scale,  since  for  a  uniformly  turbulent  path,  the  image  spreading 
due  to  turbulence  only  is  proportional  to  this  function  of  range  (c.f.  Eq. 
(24)).  Brackets  indicate  +  la  values  in  the  measured  widths.  The  solid  line 
brackets  show  standard  devTation  values  for  the  means  of  all  the  data.  Broken 
line  brackets  (1978  only)  show  standard  deviations  for  a  single  datum. 

Adjacent  numbers  give  the  number  of  frames  of  data  used  in  finding  the  means 
and  standard  deviations.  Circles  show  the  mean  spread  function  widths  to  be 
expected  from  computations  of  Cj}  and  the  resulting  in  optical  propagation 
(c.f.  Eq.  (24)).  The  zero  range  value  was  obtained  by  Fourier  transforming 
the  combined  camera  and  film  responses  of  (Eq's.  (21)  and  (22)).  The  ratio  of 
observed  to  meteorologically  predicted  line  widths  for  each  range  (excluding 
the  data  for  the  500  m  range*)  have  been  averaged  geometrically  for  each  of 
the  figures.  This  average  ratio  is  represented  on  the  figures  byX. 

b.  Factors  Concerning  Specific  Data  Groups 

Table  II  summarizes  factors  peculiar  to  the  data  groups. 


*  lie  do  this  because  data  from  the  1978,  1980,  and  1981  observation  series 
show  excessive  line  widths  for  this  near  range.  We  must  conclude  that 
factors  other  than  the  atmosphere  lead  to  this  line  broadening. 
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Atmospheric  Stability 


— 2  /  ^ 

C2  ((maters)  ) 


(at  6.7m.  elev. ) 


Wave  Height  (m. ) 


Absolute  Humidity  (gm. /kg.) 
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FIGURE  11 

VALUES  OF  METEOROLOGICAL  INPUT  VARIABLES  MEASURED  ON  THE 
YFRT-287  AND  RESULTING  COMPUTED  VALUES  OF  Cn — JUNE  6,  1980 
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FIGURE  12 

METEOROLOGICAL  DATA  TAKEN  ON  MAY  8,  1981 
AND  RESULTING  COMPUTED  VALUES  OF  Cg 


28 


RflME  E:  <  M  .  > 

13. a  THEORETICAL  (CIRCLES)  AND  MEASURED  LINE  WIDTHS  FOR 
OBSERVATION  SERIES  A  AS  A  FUNCTION  OF  TARGET  RANGE 
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Figure  13.b  THEORETICAL  (CIRCLES)  AND  MEASURED  LINE  WIDTHS  FOR 
OBSERVATION  SERIES  A  AS  A  FUNCTION  OF  TARGET  RANGE 


Figure  14. a  THEORETICAL  (CIRCLES)  AND  MEASURED  LINE  WIDTHS  FOR 
OBSERVATION  SERIES  B  AS  A  FUNCTION  OF  TARGET  RANGE 


14. b  THEORETICAL  (CIRCLES)  AND  MEASURED  LINE  WIDTHS  FOR 
OBSERVATION  SERIES  B  AS  A  FUNCTION  OF  TARGET  RANGE 


RUN  5ER  I  ETS 


gure  14. C  THEORETICAL  (CIRCLES)  AND  MEASURED  LINE  WIDTHS  FOR 
OBSERVATION  SERIES  B  (CONT’D)  AS  A  FUNCTION  OF  TARGET  RANGE 


Figure  15. a  THEORETICAL  (CIRCLES)  AND  MEASURED  LINE  WIDTHS  FOR 
OBSERVATION  SERIES  C  AS  A  FUNCTION  OF  TARGET  RANGE 


15. b  THEORETICAL  (CIRCLES)  AND  MEASURED  LINE  WIDTHS  FOR 
OBSERVATION  SERIES  C  AS  A  FUNCTION  OF  TARGET  RANGE. 


16. a  THEORETICAL  (CIRCLES)  AND  MEASURED  LINE  WIDTHS  FOR 
OBSERVATION  SERIES  F  AS  A  FUNCTION  OF  TARGET  RANGE 
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16. b  THEORETICAL  (CIRCLES)  AND  MEASURED  LINE  WIDTHS  FOR 
OBSERVATION  SERIES  F  AS  A  FUNCTION  OF  TARGET  RANGE 


Figure  17. a  THEORETICAL  (CIRCLES)  AMD  MEASURED  LINE  WIDTHS  FOR 
OBSERVATION  SERIES  G  AS  A  FUNCTION  OF  TARGET  RANGE 


Figure  17. b  THEORETICAL  (CIRCLES)  AND  MEASURED  LINE  WIDTHS  FOR 
OBSERVATION  SERIES  G  AS  A  FUNCTION  OF  TARGET  RANGE 


FIGURE  18 

THEORETICAL  (CIRCLES)  AND  MEASURED  LINE  WIDTHS  FOR 
OBSERVATION  SERIES  H  AS  A  FUNCTION  OF  TARGET  RANGE 


FIGURE  19 

THEORETICAL  (CIRCLES)  AND  MEASURED  LINE  WIDTHS  FOR 
OBSERVATION  SERIES  I  AS  A  FUNCTION  OF  TARGET  RANGE 


MEASURED  LINE  WIDTHS  FOR  OBSERVATION 
SERIES  K  AS  A  FUNCTION  OF  TARGET  RANGE 


FIGURE  22 

THEORETICAL  (CIRCLES)  AND  MEASURED  LINE  WIDTHS  FOR 
OBSERVATION  SERIES  L  AS  A  FUNCTION  OF  TARGET  RANGE 


FIGURE  24 

THEORETICAL. (CIRCLES)  AND  MEASURED  LINE  WIDTHS  FOR 
OBSERVATION  SERIES  N  AS  A  FUNCTION  OF  TARGET  RANGE 
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VII.  DISCUSSION  OF  RESULTS 

a.  Unstable  Atmosphere  with  Strong  Turbulence 

Under  these  conditions  we  note  that  mean  values  of  measured  line  widths 
differ  at  most  by  a  factor  of  two  from  values  predictable  from  bulk 
meteorological  measurements.  This  Is  in  approximate  agreement  with  the  factor 
of  two  accuracy  which  Davidson,  et  al la, 4  obtained  the  index  structure 
constant  from  bulk  meteorological  measurements.*  It  Is  significant  that  the 
geometric  mean  of  the  ratio  of  observed  to  predicted  line  widths  for  the  full 
set  of  data  Is  much  closer  than  this,  viz.,  0.91/1,  indicating  that  when  space 
and  time  variations  are  averaged  out,  the  correspondence  of  theory  with  mean 
turbulence  conditions  is  much  closer.  A  careful  look  at  the  detailed  records 
of  ambient  conditions  (Figures  6  present  30  minute  averages)  shows  that  from 
12:00  to  15:00  on  September  26,  the  wind  speed  was  highly  variable.  values 
for  data  series  B  and  C  (Figures  15  and  16)  for  the  intervals  11:45-12:50  and 
14:00-14:30  on  this  day  are  0.85  and  0.60:  giving  a  geometric  mean  of  0.71; 
whereas  the  values  for  the  other  data  series.  A,  F  and  G,  are  1.17,  1.01,  and 
1.04,  giving  a  geometric  mean  of  1.07. 

Concerning  single  line  width  observations— an  aspect  which  has  not  been 
treated  theoretically --we  note  that  the  lo  fractional  deviation  from  the  mean 
value,  averaged  over  the  set  of  ranges  (excluding  460  m. )  and  times  was  34 
percent  of  the  mean.  The  maximum  la  value  at  maximum  range  was  86  percent. 


b.  Neutral  Atmosphere 

In  this  series  of  observations  the  measured  line  width  values  are  higher 
than  the  calculated  values  and  essentially  independent  of  target  range.  This 
points  to  an  instrumental  limitation  in  the  achievable  line  width,  apparently 
about  25  to  30  microns. 

In  Figure  18,  a  correlation  of  measured  line  width  with  computed  line 
width  Is  shown  at  ranges  of  2990  and  4550  m.  At  these  times,  the  products  of 
C fj  (computed  for  a  height  of  6.67  m. )  and  ranges  were  8.0  and  2.1  x 
10-14  m.1/3.  From  this  data  one  concludes  that  if  an  ideal  camera  located 
one  meter  above  the  sea  surface  records  Images  of  a  ship  or  shore  in  a 
millisecond  exposure,  the  images  will  have  a  mean  line  spread  of  2.5  arc 
seconds**  or  less  if  the  two  primary  ambient  conditions  are  met: 

1.  Air  Temperature  -  Water  Temperature  £  1.25°C  (stable  atmosphere). 

2.  Wind  Speed  2.0  m./sec.  *  3.9  knotsT 


*  A  factor  of  two  agreement  In  line  width  can  be  shown  to  be  equivalent  to  a 
factor  of  three  agreement  in  the  values  of  cf}  . 

**  This  value  is  obtained  by  rss  subtraction  of  the  limiting  value  of  the 
line  width  produced  by  the  equipment  from  the  mean  measured  line  width 
values.  It  also  corresponds  approximately  with  the  line  width  predicted 
for  an  Ideal  camera  with  a  x  range  product  of  8.0  x  10-1* 
m.1/3  with  a  stable  atmosphere. 
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As  a  spot  check  on  the  range  dependence  of  the  mean  line  widths,  we  noted 
that  in  Figures  16b  and  17b,  where  the  temperature  and  wind  conditions  were 
stabilized,  the  (range)0*6  dependence  of  the  line  width  was  followed  with  a 
deviation  of  less  than  6  percent  over  the  range  span  920  to  2740  m.  As  a  spot 
check  on  the  camera  elevation  factor,  we  compared  two  measured  line  width 

values  at  a  range  of  2740  n.  With  Cft  values  nearly  the  same,  lowering 
the  camera  from  a  height  of  2.3  m.  (Figure  16)  to  0.8  m.  (Figure  17)  increased 
the  line  width  by  47  percent;  whereas  theoretically  this  change  (allowing  for 
different  L  values)  should  produce  a  change  of  only  40  percent.  The 
difference  ratio,  1.47/1.40  *  1.05,  is  so  close  to  1.0  that  we  conclude  that 
the  camera  height  formulation  is  correct. 

Because  of  the  indicated  correspondences  between  measured  mean  and 
theoretical  values  of  line  spread,  we  conclude  that  the  theoretical 
formulation:,  are  accurate  within  tens  of  percent  for  these  turbulence 
conditions  for  target  ranges  from  910  to  2740  m.  and  wind  speeds  from  2  to  7 
m./sec.  and  camera  elevations  from  the  surface  of  the  ocean  from  0.8  to  2.4  m. 

c.  Stable  Atmosphere  —  Weak  Turbulence 

Since  the  particularly  high  air-water  temperature  differences  shown  in 
Figure  12  in  the  time  period  09:45-11:00  do  not  lead  to  large  measured  line 
widths  in  Figure  21  ( 10:05-1 1:18),  we  conclude  that  in  this  time  interval  the 
backup  air  temperature  measuring  thermometer  was  subject  to  localized  heating 
and  thus  was  not  measuring  the  true  temperature  of  the  air  above  the  ocean. 
From  the  graph  of  wind  speed  vs-  time  we  note  that,  with  the  exception  of  a 
brief  burst  at  10:30,  the  air  speed  was  less  than  1  m./sec.  in  this  time 
i nterval . 

In  Figure  22  (11:30-12:15)  the  line  widths  recorded  appear  to  fall  in  the 
region  of  the  minimum  possible  with  the  instrumentation,  with  the  possible 
exception  of  the  30.6  micron  line  width  obtained  at  a  target  range  of  5000  m. 
at  11:30.  Generally,  the  line  width  values  computed  theoretically  appear  to 
be  unduly  high.  Again,  this  may  be  the  result  of  the  low  wind  speed,  1.3 
m./sec. 


In  Figure  23  (12:35-14:00)  the  wind  had  stabilized  at  2.5  to  3.5  m./sec., 
but  the  measured  air-water  difference  values  were  from  1/2  to  1#C  too  high  to 
permit  agreement  of  measured  and  computed  line  width  values. 


The  measured  and  computed  line  width  values  shown  on  Figure  24  (excluding 
the  datum  at  500  m. )  show  reasonable  agreement.  Although  the  theoretical 
values  are  based  on  only  one  meteorological  measurement  set,  the  14:30 
recording,  the  wind  and  temperature  had  stabilized  by  this  time.  Thus,  the 
extrapolation  of  meteorological  data  to  15:00  seems  justified. 

We  assume  that  the  results  of  Figures  23  and  24  are  based  on  reasonably 
correct  data  and  uniform  ambient  conditions.  Then  the  combined  ratio  of 
measured/computed  line  widths  for  the  figures  is  (0.61  x  1. >1/2  =  0.78. 
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VIII.  GENERAL  CONCLUSIONS 

For  strong  turbulence  with  unstable  atmospheric  equilibrium,  mean 
measured  line  widths  will  match  the  values  computed  from  bulk  meteorological 
measurements  by  the  formulations  given  In  this  monograph  to  within  a  few 
percent.  Single  line  width  measurements  can  differ  from  the  computed  values 
by  a  factor  of  up  to  two  times. 

With  neutral  or  weakly  stable  atmospheric  conditions  (Cjj  <  iq-14 
m.2/3),  rainlnum  line  width  values  varying  from  3.5  to  4  arc  sec.  were 
recorded,  rather  than  the  stable  2.5  arc  sec.  line  width  expected.  This 
seriously  limited  the  accuracy  with  which  results  could  be  related  to  computed 
line  width  values,  and  thus  the  accuracy  with  which  the  effects  of  the 
turbulence  could  be  inferred.  Two  sets  of  mean  line  widths  obtained  under 
atmospherically  stable  conditions  gave  combined  ratios  of  measured  to  computed 
line  widths  of  0.78/1.0,  but  the  turbulence  contributions  to  the  line  widths 
were  too  small  to  draw  conclusions  about  the  accuracy  of  the  theoretical 
turbulence  formulations. 
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APPENDIX  I 


COMPUTATION  OF  Cfl(z)  FROM  BULK  METEOROLOGICAL  MEASUREMENTS 


From  initial  bulk  meteorological  measurements  made  at  a  height  z'  we 
proceed  to  derive  C£(z),  where  z  is  any  height  above  the  surface.  This 
will  be  done  for  both  stable  and  unstable  atmospheres. 

The  first  variable  to  be  obtained  is  the  surface  roughness  distance  zo 
given  by 

20  =  10  exp(-0.35//CoN) , 

where  Cqn  is  obtained  from  the  wind  speed  U  at  z1  from  the  tabulation 


Con  x  1000 

1.080  U"0,15 
0.770  +  .086U 
0.870  +  .067U 
1.200  +  0.25U 


U(m/s) 

0.3  to  2.2 
2.2  to  5.0 
5.0  to  8.0 
8.0  to  25.0 


This  leads  to  values  of  the  drag  coefficient  Cdn, 
CDN(U)  =  [0.35/ln(z7zo)]2. 


(A-l ) 


We  then  calculated  a  stability  correction  factor  Co»  which  is  used  in 
setting  up  scalar  profiles  of  wind  speed,  temperature,  and  humidity.  This  is 
given  by 


,  1.62z  •  (a&*6.2x10  i a  ) 

=  _ 9 _ 


T  .  In(5xl0  z‘ )  •  CDN(U).U‘ 


(A-2 ) 


where  the  variables  are  measured  at  height  z'. 

From  the  stability  correction  factor  the  true  atmospheric  stability  at 
height  z'  is  given  implicitly  by 


K  *  z  /l*ZQ  • 


[l-i^C  )/ln(zA0)]2 
l-^(c')/ln(5xl04z) 


(A-3a) 


For  stable  atmospheric  conditions  (Co  >  0)* 
^(C  )  *  -4.7c' 


(A-3b) 


^2(C')  -  -6.5  C 
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and  for  unstable  atmospheric  conditions  (Co  <  0), 

^(C1)  =  2  ln[(l  +  x)/2]  +  ln[(l  +  x2)/2]  -  2  tan_1(x)  +  */2, 


wi  th 


x  =  (1  -  15  c' ) 1/4 
and 


(A-3c ) 


^(C)  =  2  In  [(1  +  Si  -  9C')/2] 

For  a  stable  atmosphere.  Equations  (A-3a)  and  (A-3b)  nay  be  formulated  as  a 
quadratic  function  of  C'(z')  and  solved  for  C‘ •  For  an  unstable  atmosphere 
Equations  (A-3a)  and  (A-3c)  may  be  solved  iteratively  for  C'  by  substituting 
first  C'o  and  then  successively  more  accurate  values  of  C1  in  the  right  hand 
side  of  Eq.  (A-3a)  and  calculating  the  resulting  improved  value  of  c  after 
each  substitution. 

We  now  require  values  for  the  temperature  and  humidity  scaling  parameters, 
T*  and  Q*  (gn/m3),  respectively.  They  are  given  by 

T*  =  0.47A9/[ln(5  x  104z‘)  -  ^U')]  (A-4) 


and 


Q*  =  1.3  q*  *  0.61  aq/[ln{ 5  x  10*z')  -  ^ 2f  C* )J  (A-5) 

From  this  point  on  we  refer  to  the  height  z  at  which  the  turbulence  will 
be  observed.  Then  the  stability  parameter 

C(z)  *  z/L  *  Jr  *  E'(z) '  (A-C) 


From  initial  temperature  and  humidity  difference  data  and  the  derived 
value  of  C  we  may  now  calculate  structure  functions  for  temperature  and 
humidity,  as  well  as  the  temperature-humidity  covariance  function.  These  are 
given  by 

C2  =  T2.z~2/3-f(E),  (A-7a) 

C2  *  0.6  Q2*z”2/3*f(C),  (A-7b) 

CTQ  -  0.62  T*Q*.z"2/3-m),  (A-7c) 
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where  the  height  dependence  function  f  is  given  for  a  stable  atmosphere  by 
f ( C)$  =  4.9  (1  +  2.4£2/3) 
and  for  an  unstable  atmosphere  by 
fU)u  =  4.9  (1  -  70‘2/3  . 


Finally,  by  combining  the  effects  of  temperature  and  humidity  and  relating 
these  to  variations  in  optical  index,  we  obtain  for  the  index  structure 
functi on 


Cn2  ■  (79  x  10'6P/T2)2.(C2  +  0.113  CTQ  +  0.0031  Cq),  (A-8) 

where  Pis  the  air  pressure  in  mb  and  T  is  the  absolute  air  temperature. 

For  an  optical  path  having  a  varying  height  above  the  sea  surface,  it  is 
useful  to  have  an  expression  for  giving  an  explicit  height 
dependence.  Combining  Equations  (A-6),  (A-7),  and  (A-8)  gives 

C2  =  (79  x  10*6P/T2)2-(t2  +  0.070  T*Q*  +  0.0019  Q*1  z‘2/3-f (z/L).  (A-9) 
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